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G\ ; Abstract 

0^ , New boundary bound states (BBS) are found of an integrable model with the mag- 

netic impurities located at the edges of an open Heisenberg spin chain. These bound 
states carry the real energy and are formed by three or five imaginary modes of the 
rapidities. These imaginary modes of the rapidities give the non-zero antisymmet- 
ric wave functions and the moments of the centers of the bound states are zero. 
It means that these bound states are arisen by the magnetic impurities and local- 
ized at the edges of the correlated system. The Kondo screening occurs for the 
. antiferromagnetic spin chain with the ferromagnetic impurities-electrons exchange 

I 1 interaction. 

t— > . 

Magnetic impurities in one-dimensional (ID) strongly correlated electron systems or quantum spin chains have been 
the focus of intense investigation. These strong correlated systems can be described in term of Luttinger liquid 
and the behavior of the impurities in the ID quantum systems is rather different from that in a Fermi liquid 
The availability of nonperturbative techniques allows us to have detailed pictures understanding the relevant physics 
'-^ | and revealing some very interesting phenomena such as Kondo problem jj^] and pinning of bound states in the low 
Ch ' dimensional strong correlated electron systems and quantum spin chains. Experimentally, the magnetic impurities 
implant in carbon nanotubes or quantum wires and analogical phenomena, for example, the x— ray boundary effects, 
i metal point-contact spectroscopies, etc, renew also the investigations of these problem. 

The quantum inverse scattering method (QISM) and the Bethe ansatz (BA) techniques have provided us the very 
t-H , effective tools to study the magnetic impurities in the ID quantum system. These methods have been used successfully 
to deal with the Kondo impurity in a free electron host, the magnetic impurities in spin chains and the mixed valent 
behavior of hybridization (Anderson-like, with hybridized impurity and host wave functions) |^|-^|. Recently, in 
correlated electron hosts, the properties of the magnetic impurities have been studied in a series of very interesting 
papers [^|^]. The periodic boundary conditions were imposed on the electron host and spin chains for all these cases. 
Kane and Fisher investigated a ID repulsive interacting system in the presence of a potential barrier and pointed 
out that it corresponds to a chain disconnected at the barrier site at low energy scales fliotl . This can be effectively 
0^ ■ described by the open boundary conditions which are well studied by the boundary conformal field theory [ pd| and 
the BA methods p^JT^] . Zvyagin get that the low-energy magnetic behaviors of an impurity in a chain with periodic 
boundary conditions and with open boundary conditions coincide up to mesoscopic corrections of order of L^ 1 with L 
{~j being the length of the system || . Now we know that several methods have been used for introducing the impurities 
into the integrable models of correlated electrons and quantum spin chains with the open boundary conditions. The 
first one for the construction of the impurity models is dependent on the main idea that the spectral parameters in the 
scattering matrix of the model rely on the difference of the particles' rapidities. This method has been firstly used to 
the impurity model of the ID quantum lattice gases with the periodic conditions by Eckle, Punnoose and Rbmer M, 
where the integrable condition continues to be satisfied under an arbitrary local shift of the related parameter (See, 
for example pj). The second method depends on that the scattering matrix of the bulk of the t — J model relies on 
the tangent ( or cotangent ) functions of the half moments of the electrons [fl5l| . The very interesting thing is that 
the corresponding Hamiltonian with the magnetic impurities has a simple and compact form Jig] , which is different 
from the one in the open Hubbard impurity model |17|] . The boundary scattering matrix between the impurity and 
the electron can be factorized as two terms: one is similar as the R matrix and another is similar as the inverse of the 
R matrix with the inverse spectral parameter, justly as rightly obtained by Zvyagin and Johannesson in their very 
interesting work |l8| . They revealed the existence of a hidden Kondo effect driven by forward electron scattering off 
the impurity related to this property of the boundary scattering matrix. 

As is well known, the bound states can be formed for the strongly correlated electron system and quantum spin 
chains within the charge or spin sectors and they are very important to decide the thermodynamics of the system 
and the its low temperature properties. Despite the success of the QISM and BA approach to the investigation of the 
magnetic impurities in correlated system, it is not very clear that the impurities contribute to the bound states of the 
strongly correlated system with the open boundary conditions. In this letter, we discuss this problem in details. 
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In general a bound state can be formed by several complex rapidities of the particles, such as charges or spins. The 
total energy of these complex modes and the total moments of the ones should be real. Under the open boundary 
conditions, the rapidities Uj of the charges or spins should satisfy also that Uj ^ ±u/ when j ^ I. Otherwise, the wave 
function is zero, which means that this situation is forbidden. When the complex rapidities form a bound state, the 
above three conditions should be satisfied for the correlated system under the open boundary conditions. For clarity 
we here focus on the case of the quantum Heisenberg model. When the two impurities with the arbitrary spins are 
coupled to this open quantum spin chain, the Hamiltonian of the system can be written as 

H = — &j ' + Jl <?i • Sl + Jr&n ■ Sr, (1) 

4=1 

where dj are the Pauli matrices; Sl,r are the impurity moments with the arbitrary spins Sl,r- The site number of 
the bulk is N . Jl,r are two arbitrary real constants which describe the coupling between the bulk and the impurities 
and can be parameterized as 

J l,R = J — 2 — (2) 



(Sl,. 



2 J 'L.I! 



with the arbitrary constants cl.r. This Hamiltonian can be diagonalized by using the standard Bethe ansatz scheme. 
The eigenvalue of the energy of this impurity system is 



M 



E(Xi,X 2 , ■ ■ ■ , A M ) - ^ j i + ^ JiSi-\ — -. 



j=l 1 4 l=L,R 



with the following Bethe ansatz equations: 




Xj +i(Si+ rci) 
Xj -i(Si + rci) 



- n w 

1 = 1(1^3) ''=±1 3 

Now we study the boundary bound states in detail for the above Heisenberg impurity model. From definition (||) we 
know that —Cl,r is equivalent to cl,r because they give the same Hamiltonian (1). So, without losing generality, we 
restrict the parameters cl.r to non-negative values in the following discussion. 

When the coupling in the bulk is antiferromagnetic ( J > 0) and the coupling between the bulk and the impurity is 
ferromagnetic, the system has the following boundary bound states (BBS): 

A3,i = i(Sl,r - cl,r), (4) 
A3, 2 = —^(Sl,r - C L .R ■- 1), (5) 

2 

A3.3 = —^(Sl,r - c l ,r + 1), (6) 

when the parameters cl.r in the regime Sl.r. + 1/2 < Cl,r < Sl,r + 1. These two BBS, one BBS for each ends ( L 
and R) of the chain, are formed by the three imaginary modes of A, respectively. They carry the energy 



p(3) 12 J [3(S L , R - c L , R f - 2] 

L R [±(Sl,r - c L , R y - 1] [{S L , R - c L , R Y - 4] • 



The moments of the centers of the BBS are Ylj=i ^3,j = 0- ^ means that this kind of bound states is localized at 
the edges of the system. Of course, the energy and the moments of the centers of the BBS are all real and the spin 
rapidities satisfy that A3_i =/= ±A3^ =/= A3. 3 ^ iAs^ which ensure that the antisymmetric wave functions of the system 
are not zero. Therefore, these BBS satisfy all the physical demand. By making the transformations A3.J — > — A3.J 
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(j = 1,2,3) in the relations (f|j6|), the BBS with the three imaginary modes can be obtained also for the Heisenberg 
impurity model with Sl,r + 1/2 < cl.r < Sl,r + 1. They carry the energy also as the form (Q). 

When the coupling in the bulk is antiferromagnetic (J > 0) and the coupling between the bulk and the impurity 
falls also into the antiferromagnetic regime, the above BBS with the three imaginary modes are formed under the 
condition I/3 + Sl.r < Cl,r < 1/2 + Sl.r- The corresponding imaginary modes with the transformations X^ j — > —Xaj 
(j = 1, 2, 3) are also the BBS in this case. The energy expression (Q) do not change the forms and the moments of the 
centers of the bound states are also zero. By taking the logarithm of the Bethe ansatz equation (^)and introducing the 
distribution functions of the spin rapidities, we can get the integral equations of the impurity model for the ground 
state. Then we have that the self-magnetization of the ground state is Sl + Sr ~ 1 for the two up impurity spins, or 
1 — Sl — Sr for the two down impurity spins, or ± (Sl — Sr) for the one up and one down impurity spins. The similar 
procedure gives that the self-magnetization of the ground state is also ± (Sl — Sr) when the coupling in the bulk is 
antiferromagnetic but the coupling between the bulk and the impurity is ferromagnetic. When the parameters Cl,r, 
which describe the coupling between the impurities and the bulk, satisfy that Sl.r < cl,r < Sl.r. + 1/3, there are 
the BBS formed by the three imaginary modes as relations and the ones by the transformations X^j — > — \$ t j 

(j = 1, 2, 3). They carry the energy (0) and have the zero moments of the centers of the BBS. In this case, the coupling 
between the bulk and the impurities are in the antiferromagnetic regime if the exchange interaction in the bulk is 
antiferromagnetic, too. Furthermore, in this case (Sl,r < cl.r < Sl,r + 1/3), there are two bound states (one for 
each ends of the chain) formed by the following five imaginary modes of A, respectively, 

As,i = i(Sh,R — cl,r), 
= i(SL,R - cl,r + 1), 

^5,3 = i(S L ,R - C L ,R - 1), 

1 

^5A = —-^(3Sl,r — 3cl,r — 1), 

1 

Xz,5 = --^(^Sl.r - 3c l ,r + 1). 

They carry the energy 

(5 ) 20 J [7(S L . R - c L m) 2 - 6] 

L ' R [9(Sl,r ~ c l ,rY ~ 4] [4(S l ,r - c l ,r) 2 - 9 

The moments of the centers of above bound states formed by the five imaginary modes of A are 2 ■■■ 5 = 0- 

So they localize at the two edges of the Heisenberg spin chain and A5J ^ ±A 5 j (j, I = 1, 2, • • • , 5) which ensure 

that the system has a non-zero antisymmetric wave function. Similarly as the case of the BBS of the three imaginary 
modes, the transformations A 5 j — > — As j (j — 1,2, • • ■ ,5) give also the boundary bound states of the system and 
they do not change the expression of the energy and the moments of the centers of the bound states are also zero. 
By the use of the method mentioned above, we get that the self-magnetization of the model is same as the situation 
1/3 + Sl.r < Cl,r < 1/2 + Sl.r- The above BBS satisfy the three conditions of the complex modes. 

In the following part, we describe simply the self-magnetization and the BBS of the ferromagnetic Heisenberg 
impurity model. By solving the Bethe ansatz equations in the thermodynamic limit, we get that the self-magnetization 
of the system is Sl+Sr+N/2 when the coupling between the impurities and the bulk is also ferromagnetic. Otherwise, 
the self-magnetization is — Sl — Sr + N/2 for the antiferromagnetic exchange interaction between the impurities 
and the bulk. When one impurity has the ferromagnetic interaction with the bulk and another impurity has the 
antiferromagnetic interaction with the bulk, the self-magnetization of the system has the form ± (Sl — Sr) + N/2. 
When Sl.r. + 1/2 < cl r < Sl,r + 1, the system has the BBS formed by the three imaginary modes (Q-^|)and the 
BBS carry the energy (j?]). By making the transformations X^j — > — X3J (j — 1, 2, 3), the corresponding bound states 
are also the BBS of the system. They satisfy all of the three conditions and the coupling between the bulk and 
the impurities is antiferromagnetic. When the coupling between the bulk and the impurities is ferromagnetic, the 
system has the BBS which can be formed by the three imaginary modes fljq) or the five imaginary modes 
with Sl.r < Cl,r < Sl.r. + 1/3. They carry the energies as expressions and ([l3]), respectively. Of course, the 
corresponding imaginary modes with the transformations Xaj = —Xa,j (j = 1) 2, 3) and As,; = — Ag; (I = 1, 2, • • • , 5) 
form also the BBS of the impurity model. When 1/3 + Sl < Cl < Sl + 1/2, we have the BBS or the inverse of 
rapidities and the exchange interaction between the bulk and the impurities is ferromagnetic. Notice that the system 
may have the other form of the impurity bound states in the above restricted range of the impurity couplings Jl,r- 
It is also an open problem to find out the impurity bound states for other values of the impurity couplings for the 
strongly correlated system. 
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(9) 
(10) 

(11) 

(12) 



The above discussion shows that the ferromagnetic and antiferromagnetic Heisenberg spin chains with the magnetic 
impurities have always the BBS arisen by the impurities when the strengths of the interactions between the bulk and 
the impurities are proper. These bound states carry the real energy and the moments of the centers of the BBS are 
zero. They satisfy all of the three conditions of the imaginary modes. For the ferromagnetic Heisenberg model, the 
system has the BBS formed by the three imaginary modes and the five imaginary modes of the spin rapidities when 
the coupling between the bulk and the impurities is ferromagnetic. When the coupling between the bulk and the 
impurities is antiferromagnetic, the system has only the BBS contained the three imaginary modes of the rapidities. 
For the antiferromagnetic Heisenberg model, the system has the BBS formed by the three imaginary modes and the 
five imaginary modes of the spin rapidities when the coupling between the bulk and the impurities is antiferromagnetic. 
When the coupling between the bulk and the impurities is ferromagnetic, the system has only the BBS contained the 
three imaginary modes of the rapidities. And the Kondo screening -as predicted by Furusaki and Nagaosa exists 
for the antiferromagnetic Heisenberg model with the ferromagnetic coupling between the bulk and the impurities. At 
zero impurity couplings (Jl,r — 0), the system has impurity bound states, although trivial, which just correspond 
to the 2Sl,r + 1 spin states of the impurity. By turning on the boundary coupling, the BBS can be formed by the 
imaginary modes of the rapidities and the number of the BBS might change. The BBS can affect the ground state 
of the whole system when the boundary coupling is strong enough, which is under investigation. New properties of 
the specific heat, excited state, dressed energy, et al can be introduced due to the impurity couplings. We point out 
that the similar BBS with three or five imaginary modes contributed by the magnetic impurities can be found also 
for the strongly correlated electron systems such as the Hubbard model and the t — J model with the open boundary 
conditions in the charge sectors. The BBS carry the energy and satisfy the above three conditions of the imaginary 
modes. Finally, the way the string affects the distribution of the rapidities is also the interesting subject for the 
further investigation. 

To conclude, we have found that the BBS arisen by the magnetic impurities in a correlated host under the open 
boundary conditions. These BBS are formed by three or five imaginary modes such as charges or spins. The imaginary 
modes of the bound states due to the magnetic impurities satisfy that: (i) the total energy is real; (ii) the total moment 
of the imaginary modes is real (zero); (iii) the absolute values of the imaginary modes are different. These BBS carry 
the energy and pin at the edges of the system. And the Kondo screening exists for the antiferromagnetic Heisenberg 
model with the ferromagnetic coupling between the bulk and the impurities. 
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